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Abstract. 

The timing and duration of the reionization epoch is crucial to the emergence and evolution of structure in the universe. The 
relative roles that star-forming galaxies, active galactic nuclei and quasars play in contributing to the metagalactic ionizing 
background across cosmic time remains uncertain. Deep quasar counts provide insights into their role, but the potentially 
crucial contribution from star-formation is highly uncertain due to our poor understanding of the processes that allow ionizing 
radiation to escape into the intergalactic medium (IGM). The fraction of ionizing photons that escape from star-forming 
galaxies is a fundamental free parameter used in models to "fine-tune" the timing and duration of the reionization epoch that 
occurred somewhere between 13.4 and 12.7 Gyrs ago (redshifts between 12 > z > 6). However, direct observation of Lyman 
continuum (LyC) photons emitted below the rest frame H I ionization edge at 912 A is increasingly improbable at redshifts 
z > 3, due to the steady increase of intervening Lyman limit systems towards high z. 

Thus UV and U-band optical bandpasses provide the only hope for direct, up close and in depth, observations of the types 
of environment that favor LyC escape. By quantifying the evolution over the past 11 billion years (z <3) of the relationships 
between LyC escape and local and global parameters such as: metallicity, gas fraction, dust content, star formation history, 
mass, luminosity, redshift, over-density and quasar proximity, we can provide definitive information on the LyC escape fraction 
that is so crucial to answering the question of, how did the universe come to be ionized? Here we provide estimates of the 
ionizing continuum flux emitted by "characteristic" (L* v ) star-forming galaxies as a function of look back time and escape 
fraction, finding that at z = 1 (7.6 Gyrs ago) V uv galaxies with an escape fraction of 1% have a flux of 10 -19 ergs cm -2 s _1 
A- 1 . 
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THE ESSENTIAL ROLE OF LYC ESCAPE 

Some 0.3 Myr after the Big Bang, the adiabatic expansion of the universe caused the primordial plasma of protons 
and electrons to cool, creating a neutral gas. Recent observations show that most of the universe has since been 
reionized and provide constraints to the duration of this process. Sloan Digital Sky Survey spectra of luminous high- 
redshift quasars have black H I Gunn-Peterson troughs, indicating a mean H I fraction of > 10~ 3 at z > 6.4 when the 
age of universe was w 1 Gyr [1]. Evidence that reionzation started even earlier is provided by the polarization of the 




FIGURE 1. The results of a Monte Carlo [10] showing the 
effects of intergalactic absorption by Lyman limit systems 
on the transmission of LyC photons. The graph depicts the 
cumulative probability of having a line-of-sight transmission 
greater than that shown on the axis for LyC photons emitted 
between (880 - 910 A). The light blue, magenta, dark blue, 
green and red are contours for the redshifts z = 1, 2, 3, 4 
and 5 respectively. At z = 4, the probability of having a 
transmission of greater than 30% is 0.2 
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microwave background on large angular scales seen by the Wilkinson Microwave Anistropy Probe, which is consistent 
with an ionization fraction ~ unity at z ~ 11 when the universe was w 365 Myr old [2]. 

Complete reionization occurs when the rate of ionizing photons emitted within a recombination time exceeds the 
number of neutral hydrogen atoms. The duration of the reionization epoch depends on the initial mass function (IMF) 
of the first ionizing sources, their intrinsic photoionization rate (g), the baryon clumping factor (C = < p 2 > / <p> 2 ), 
and the fraction of ionizing photons that somehow manage to escape into the intergalactic medium (IGM) [3]. Of these 
parameters the LyC escape fraction (f e ) is the least constrained [1]. Its (often arbitrary) choice can alter conclusions 
regarding the nature and duty cycle of the sources thought to be responsible for initiating and sustaining reionization 



LyC escape plays an essential role in the formation of structure. The escape fraction of ionizing photons from 
galaxies is the single greatest uncertainty in estimating the intensity of the metagalactic ionizing background (MIB) 
over time [6]. The MIB controls the ionization state of the IGM at all epochs and may be responsible for hiding a 
non-trivial fraction of the baryons in the universe [c.f. 7, 8]. Ionizing radiation produced by star-forming galaxies is 
ultimately related to the rate of metal production by stellar nucleosynthesis [9] . The MIB intensity is a gauge of the 
feedback into the IGM of chemicals, mechanical energy and radiation by supernovae and stellar winds. 

Ionizing radiation regulates the collapse of baryons on local and global scales [11]. Photoelectrons provide positive 
feedback for star formation by promoting the formation of H~, which in turn catalyzes the production of H2; a crucial 
coolant for collapse at high-z. Negative feedback occurs when photoionization heating raises the temperature and 
inhibits stellar collapse by increasing the Jeans mass. Photodissociation of H2 is another form of negative radiative 
feedback mediated by both Lyman- Werner photons in the 912 - 1120 A bandpass and LyC photons [12]. Whether 
ionizing radiation has a positive or negative effect on a collapsing body depends on the gas density, the strength of the 
radiation field, the source lifetime and the escape fraction of LyC photons[ll, 13]. 

Quantification of the LyC escape fraction is at the frontier of reionization physics. The high opacity of even small 
column densities of H I to ionizing radiation makes the sources very faint at all epochs, but especially at redshifts z > 
3 [3]. Figure 1 from Inoue and Iwata [10], shows the likelihood of detecting LyC escape from star- forming galaxies 
becomes increasingly improbable above z > 3, due to a progressive increase with redshift in the number density of 
intervening Lyman limit and Lya forest systems. Detections above z > 4, while not ruled out, will be extremely rare. 
This favors UV and U-band optical observations in efforts to directly identify and spatially resolve those physical 
environments that allow LyC to escape. By examining at low redshift the relationship between LyC escape and the 
local and global parameters of metallicity, gas fraction, dust content, star formation history, mass, luminosity, redshift, 
over-density and quasar proximity, we seek to understand how the universe came to be ionized. 



The fundamental question is, how did the universe come to be reionized and how long did it take? Current thinking 
posits that LyC escape from the smallest galaxies powers reionization at z ~ 6, since quasars are too few in number 



[4, 5]. 



SOURCE(S) OF REIONIZATION AND THE MIB 
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FIGURE 2. Left - LFs for quasars and galaxies at z = 6. Right - Reionizing photon production rate for quasars and galaxies at 
a redshift of z = 6 and for galaxies alone at z = 7. The production rate required to maintain a fully ionized universe for clumping 
factors of 20, 30 and 40 are indicated at the top of each graph. Higher clumping factors require more photon production. See [15] 
for details. 



to sustain reionization [3, 14, 15]. However, this conclusion depends on an extrapolation of the faint end slope of the 
galaxy luminosity function (LF, -1.6 < a < -2.0), the faint end luminosity cutoff, the clumping factor (3 < C < 45) 
and escape fraction (f e ~ 0.1 - 0.8). 

The sensitivity of this conclusion to the faint end slope and the role played by the clumping factor and f e is illustrated 
in Figure 2 taken from Yan and Windhorst [15]. On the left, the LFs for quasars and galaxies are displayed for a redshift 
of z = 6. Two extreme faint end slopes are shown for the galaxies (-1.6, -2.0) and for the quasars (-1.6, -2.6). On the 
right, two panels show the cumulative reionizing photon production rate for quasars and galaxies at a redshift of z = 
6 and for galaxies alone at z = 7. Horizontal lines drawn at the top of each panel mark the critical production rate 
required to keep the universe fully ionized for clumping factors of 20, 30, 45 (higher clumping factors require more 
photon production to overcoming clump self-shielding) and assuming f e = 0.1. 

The figures show that at z = 6 the faintest galaxies dominate the LyC production and are more likely than quasars 
to maintain the universe in a fully ionized state. The case becomes less certain at z = 7 where it has been found that 
maintaining reionization requires either a "top heavy" IMF or escape fractions 0.3 < f e < 0.8, assuming 20 < C < 45 
[16, 17], although recent work by Finkelstein et al. [18] suggests 0.1 < f e < 0.5 and 3 < C < 5. It may be that there 
are not enough star-forming galaxies early on to initiate reionization [19] and that mini-quasars might be involved 
[20]. There are also indications that the initiation of reionization above z — 7 may require a "hard" spectral energy 
distribution (SED) more characteristic of quasars [19]. 

Black holes reside in the nuclei of most if not all quiescent galaxies [21, 22], so it is simplistic to characterize 
reionization as a process caused by either quasars with f e = 1 or star- forming galaxies with f e < 1 . Some fraction of 
quasars exhibit a break at the Lyman edge likely due to obscuration by host galaxies [23, 24]. The central engines of 
active galactic nuclei (AGN) have intermittent duty cycles, so the effects of previous AGN activity within an apparently 
dormant galactic environment may reduce, for a time, the local H I density aiding LyC escape. A quasar in close 
proximity to star-forming galaxies could produce a similar effect. Such considerations become increasingly important 
as the MIB budget moves from one dominated by star- forming galaxies near z ~ 6 to one with a quasar contribution 
that peaks at z ~ 2 [25]. Wide-field observation of f e from objects in extended cluster environments can be used to 
map out the relative contribution of galaxies, AGN and quasars to the MIB in the local universe and provide a means 
to assess its spatial uniformity in the redshift range z < 2.3, which cannot be probed by the ratio of He II/H I Lya forest 
lines [24]. 



LYC AND LYa ESCAPE: ENVIRONMENTS, ANALOGS AND PROXY PROSPECTS 

Reionization appears to require LyC leakage from galaxies with f e ~ 0.1, but how LyC and Lya escape from galaxies 
is somewhat mysterious. Most star-forming galaxies have mean H I columns greater than damped Lya systems 
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FIGURE 4. Surface densities as a function of 
observer's frame apparent magnitude for galaxy 
populations with redshifts between - 0.2, 0.2 
- 0.4, 0.4 - 0.6, 0.6 - 0.8, 0.8 - 1.2, 1.8 - 2.3, 
2.4 - 3.4, estimated following Arnouts[35]. There 
are 100s - 1 0,000' s of galaxies per square degree 
per magnitude with 24 > m i5oo(i+z) > ^0 f° r ea °h 
redshift interval. 
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(DLA).The optical depth at the Lyman edge for DLAs is x x<9l2 > N m 6.3 x 1(T 18 (k/9\2f = 1260 (X/912) 3 , while 
at the line core of Lya the optical depth is, x^a = Nm 6.3 x 10~ 14 = 1.26 x 10 7 (for Vd op = 12 km s _1 ). Escape from 
such large mean optical depths requires that the interstellar medium (ISM) be highly inhomogeneous. The escape of 
Lya (f a ) is aided by velocity gradients and the presence of multi-phase media [26, 27, 28, 29]. Similarly, LyC escape 
is thought to result from galaxy porosity, low neutral density, high ionization voids, chimneys created by supernovae 
or the integrated winds from stellar clusters [30], 

Exploring the possibility of a proxy relationship between f e and f a will be extremely important in the coming 
decade as JWST seeks to identify the source(s) responsible for initiating and sustaining reionization. The brightness of 
Lya emission from the first objects is expected to be much easier to detect than their rest frame UV continuum. 
Consequently, JWST will probably have to rely upon observations of Lya escape as a proxy for LyC escape. 
Unfortunately there is no guarantee that such a proxy relationship exists, because escaping Lya photons are created 
by recombining electrons freed by the LyC photons that do not escape [31] ([Lya w (2/3)Q/ a (l — f e )exp(— t)]). It is 
essential to test the proxy hypothesis at z < 3 by obtaining simultaneous observations of LyC and Lya. 

LYC ESCAPE DETECTIONS AND THE ADVANTAGES OF SPECTROSCOPY 

Siana et al. [32] have provided the most current summary of detection efforts to date. In short they have returned mixed, 
but tantalizing results, that hint at a trend for f e falling towards low-z. Spectroscopic observations hold an advantage 
over broad band imaging by providing the means to quantify ISM and IGM attenuation by H I using Lyman series 
absorption and investigating Lya escape processes. Moreover, they provide the means to detect spectroscopically, 
absorption features from species with wavelengths shortward of the Lyman edge, as was reported by Bogosavljevic 
and Steidel [33] who found evidence for the O I X 877 in a stack of 13 spectra of z = 3 galaxies. A quantitative 
assessment of the evolving contribution of galaxies to the MIB will likely require spectroscopic surveys over wide 
angular fields to acquire the large number of observations needed for establishing LyC luminosity function [34]. 

LYC DETECTION REQUIREMENTS 

GALEX has shown there are thousands of far-UV emitting galaxies per square degree down to its limiting magnitude 
mpuv ~ 25.McCandliss et al. [36] have suggested a wide field spectroscopy survey as an efficient way to search for 
LyC and Lya leakage. Instrument requirements can be derived from the sensitivity required to detect LyC escaping 
star-forming galaxies over the redshift interval 0.02 < z < 3. 1 



1 The lower limit is set by the need to work at redshifts high enough to escape the H I "shadow" of the Milky -Way. 
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FIGURE 5. The purple asterisks show the characteristic apparent LyC magnitudes (ab) w* 90 o(i+ z ) as a function of look back 
time, and in redshift and wavelength space, for different escape fractions. Contours of constant flux units are overplotted as green 
dashes marked in FEFU fractions; the background limit for FUSE. See [36] for details. 



To estimate the sensitivity requirement we use the surface density of UV rest frame emitting galaxies as a function 
of apparent magnitude (observer frame flux) and redshift, Figure 4, where the LFs from Arnouts et al. [35], for the 
redshift intervals indicated in the caption, are shown. Each asterisk marks the characteristic magnitude of the LF, 
appearing in the Schechter function. 

We convert the 1500 A characteristic apparent magnitudes to LyC magnitudes using, 

m 900(l+z) = m 1500(l+z) + ^ m 900° + Sm e , (1) 

with bm e = —2.5logf e and S^qq = 2.51og(/i5oo//90o)- Starburst99 models[37] for continuous star-formation, 
assuming solar metallicity, a Salpeter IMF and an upper mass cutoff of 100 M , yield /1500//900 ~ 2. This ratio 
is insensitive to age with 1.5 < /1500//900 ^ 3 for ages 10 - 900 Gyr. The apparent magnitude in the LyC, as a function 
of z for f e = 0.01, 0.02, 0.04, 0.08, 0.16, 0.132, 0.64 and 1 are displayed in Figure 5 as a series of purple connected 
asterisks. Dashed green lines give the conversion from magnitude to flux units. We find that L* v galaxies with f e = l% 
have LyC fluxes < 10 -19 ergs cm -2 s _1 A -1 at look back times between 7.6 and 11.2 Gyrs (1 < z < 3). 

CENTRAL QUESTIONS, PRIORITIZATION AND ENABLING TECHNOLOGIES 

Understanding the mechanisms of reionization hinges on understanding how f e changes as a function of luminosity 
and redshift. The answer will be important regardless of the outcome. If star-forming galaxies are found with f e > 0.1 
then they become plausible sources of reionization. If not, then new physics may be required to explain reionization 
[38]. 

Four central questions of fundamental importance to the field can be addressed in the coming decade: 

1. What are the relative contributions of quasars, AGN and star-forming galaxies to the MIB at z < 3. 

2. What is the relationship between f e and the local and global parameters of metallicity, gas fraction, dust content, 
star formation history, mass, luminosity, redshift, over-density and quasar proximity. 

3. Do low-z analogs exist of the faint high-z galaxies thought to be responsible for reionization? 

4. Can the escape fraction of Lya photons serve as a proxy for f e l 
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The answer to the last question is critical to the JWST key project seeking the source(s) of reionization. 

These questions can be most efficiently addressed using ultra- sensitive wide field spectroscopic surveys from space 
to probe the ionizing characteristics of galaxies as old as 1 1 Gyrs. Key mission-enabling technologies that will support 
the development of 0.°5 field-of-view, multi-object UV spectroscopy include microshutter arrays [39], high efficiency 
aberration corrected dual-order gratings [40], large format detectors [41] and GaN photocathodes [42]. Future high- 
grasp UV telescopes sensitive enough to detect LyC leak will also easily detect the cosmic web of low surface 
brightness Lya emission, providing an unambiguous beacon for emerging complexity [43, 44]. 
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